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Abstract 

, Recently, the MILAGRO collaboration reported on the detection of a diffuse multi- 

TeV emission from a region of the Galactic disk close to the inner Galaxy. The 
emission is in excess of what is predicted by conventional models for cosmic ray 
propagation, which are tuned to reproduce the spectrum of cosmic rays observed 
locally. By assuming that the excess detected by MILAGRO is of hadronic origin 
and that it is representative for the whole inner Galactic region, we estimate the 
. expected diffuse flux of neutrinos from a region of the Galactic disk with coordinates 

—40° < I < 40°. Our estimate has to be considered as the maximal expected 
neutrino flux compatible with all the available gamma ray data, since any leptonic 
contribution to the observed gamma-ray emission would lower the neutrino flux. The 
diffuse flux of neutrinos, if close to the maximum allowed level, may be detected 
by a km 3 -scale detector located in the northern hemisphere. A detection would 
unambiguously reveal the hadronic origin of the diffuse gamma-ray emission. 
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1 Introduction 



Recently, the MILAGRO collaboration reported on the detection of multi-TeV 
diffuse gamma-ray emission from two regions in the Galactic plane: the Cygnus 
region (located at galactic longitude 65° < I < 85°) pQ and the portion of the 
inner Galaxy visible from the location of MILAGRO (30° < I < 65°) 0. When 
placed alongside the earlier lower energy measurements by EGRET [3J, such 
detections provide an interesting insight into both the spatial distributions of 
cosmic rays (CR) in the Galaxy and their energy spectrum. 

The diffuse GeV emission detected by EGRET along the galactic plane is be- 
lieved to be the result of hadronic interactions between CRs and interstellar 
matter. A strong point in favor of this hypothesis is the fact that the diffuse 
GeV emission correlates with the spatial distribution of the interstellar gas 
that constitutes the target for CR interactions [3J . The same argument might 
be extended also to the diffuse multi-TeV emission detected by MILAGRO, 
since both the Cygnus region and the inner Galaxy are characterised by an 
enhancement in the gas density. In the case of the Cygnus region, the mor- 
phology of the multi-TeV emission correlates with the CO emission, which 
traces the gas density. This suggests again a hadronic origin also for the very 
high energy diffuse emission pQ. Similarly, for the region close to the inner 
Galaxy, the rather narrow latitude profile of the multi TeV-emission, concen- 
trated around the dense Galactic disk, might suggest a hadronic interpretation 
of the gamma-ray data, though also a leptonic one seems feasible [TU1[TT1I2"] . 

The expected hadronic gamma-ray emissivity from the Galactic plane can be 
calculated by solving the transport equation describing the propagation of 
CRs in the Galaxy and using the measured distributions of the gas density 
in the Galactic disk (see [I] for a review). Conventional propagation models 
are tuned in order to match the CR spectrum observed locally [5H6fT] . These 
models reproduce quite fairly the observed level of the GeV emission (with 
the notable exception of the so called GeV-excess) [8] but fail to reproduce 
the MILAGRO observations by about one order of magnitude in the case 
of the Cygnus region pQ and a factor of ~ 5 for the inner Galaxy [2J. One 
possible solution of the problem is to assume that the average CR spectrum 
in the Galaxy is harder than the one we measure at Earth [9fT0] . Another 
possibility is to assume a hard spectrum for the CR electrons [TUfTI] or that 
the typical CR electron intensity in the Galaxy is significantly higher than 
that measured locally [5]. This would make the inverse Compton contribution 
to the total gamma-ray emission dominate over the hadronic one at the multi- 
TeV energies relevant here. Remarkably, both these scenarios (hadronic and 
leptonic) were proposed to solve the problem of the GeV-excess and can now 
be generalised to investigate their predictive power at multi- TeV energies. 
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Thus, the origin of the multi-TeV emission from the central Galactic plane re- 
gion remains unclear. Either leptonic or hadronic processes are believed to be 
the dominant contributors to the detected gamma-ray flux, with the possibil- 
ity of a transition from a dominant hadronic contribution to a leptonic one in 
the GeV-TeV energy range [8]. In the case of the hadronic scenario, one would 
expect, as it is observed, a narrow extension in latitude of both the EGRET 
(see Fig. 3 from [3]) and MILAGRO emission [2J. In this case the extension of 
both GeV and multi-TeV emission is determined only by the gas distribution 
in the Galactic disk, which constitutes the target for proton-proton interac- 
tions. In the leptonic scenario, the multi-TeV emission is produced via inverse 
Compton scattering by « 100 TeV electrons which can propagate over a dis- 
tance of ~ 100 pc before being cooled by synchrotron and inverse Compton 
losses [10]. Thus, also in this case a narrow latitude distribution of gamma rays 
is expected, if the sources of cosmic ray electrons are concentrated around the 
galactic plane. In this case, however, the extension of the gamma-ray emission 
would be energy dependent, since higher energy electrons are cooled faster 
and can propagate shorter distances [10] . Thus, an accurate comparison of the 
extension in latitude of the emission detected at GeV and multi-TeV energies 
is of crucial importance, and future observations with improved angular reso- 
lution by GLAST and HAWK might help in discriminating between the two 
scenarios. 

Since the production of gamma rays is accompanied by a corresponding neu- 
trino flux only in the hadronic scenario, neutrino telescopes provide a unique 
opportunity to disentangle the origin of the mufti- TeV radiation. In this paper, 
we investigate the capability of km 3 -scale neutrino detectors such as IceCube 
[12] and KM3NeT [13J, to detect the diffuse flux of neutrinos from the in- 
ner Galaxy and possibly reveal the origin of the diffuse emission detected by 
MILAGRO. 

In the past, several estimates of the diffuse neutrino flux from the Galactic 
plane region have been published. The simplest of these approaches assumes 
that the CR distribution is constant across the Galactic disk [bf|15lll6j . whilst 
more sophisticated calculations assume a distribution of CR sources in the disk 
and obtain a steady state CR distribution by solving the diffusion equation 
[6 7J. All of these approaches, however, result in relatively small neutrino fluxes 
and are below that required to be detectable by a km 3 -size neutrino detector 
such as IceCube or KM3NeT, within one year of data acquisition. In these 
calculations, a CR spectral index of ~ 2.7 is generally assumed, or derived 
by solving the CR transport equation, since this matches the observed CR 
spectrum measured at Earth. Such a steep spectrum, normalized to reproduce 
the diffuse emission observed at GeV energies by EGRET, naturally results 
in a conservative estimate of the diffuse neutrino flux at multi-TeV energies 
where km 3 neutrino telescopes are optimally sensitive. However, no a priori 
reason exists that such a CR spectrum is applicable throughout the Galaxy. 
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Indications exist for a departure from the observed local CR spectral index 
of ~ E~ 2 ' 7 in the diffuse gamma-ray emission observed by H.E.S.S. from the 
Galactic centre ridge [T7| and the GeV excess in the diffuse flux from the 
Galaxy observed in EGRET data, which might be due to a CR spectrum 
harder than the one observed at Earth [91T10] . 

In this work we adopt a phenomenological approach. We assume that the 
production of the multi-TeV diffuse flux observed by MILAGRO is entirely due 
to hadronic interactions. This model requires the assumption of a CR spectrum 
in the inner part of the Galaxy that is different (harder) than the canonical 
2.7 spectral index of the CRs observed at Earth. Consequently we obtain a 
higher value for the multi-TeV neutrino flux from the same Galactic plane 
region compared to previous canonical estimates. Any leptonic contribution 
to the multi-TeV diffuse emission will lower the actual diffuse neutrino flux. 
Our calculations therefore provide an estimate of the maximal neutrino flux 
that can be expected, under the constraint of remaining compatible with all 
the presently available gamma-ray data. Since the MILAGRO detection of the 
diffuse multi-TeV gamma-rays from the inner Galaxy only covers the angular 
range 30° < I < 65° 0, to apply this MILAGRO data point to the whole 
inner Galaxy with coordinates —40° < I < 40° it was assumed that the TeV 
flux detected by MILAGRO scales in a similar way to the GeV flux detected 
by EGRET from the two regions. 

It is worth noticing that MILAGRO is operating at multi-TeV energies which 
is also the energy range where km 3 -scale neutrino telescopes reach their best 
performances (see [T8lfT9] and Sec. 13.11) . This implies that the MILAGRO ob- 
servations are the most relevant gamma-ray observations to constrain the ex- 
pected diffuse neutrino flux, since gamma rays and neutrinos are produced 
roughly at the same energy during CR hadronic interactions with the inter- 
stellar medium. 

In Sec. [2]we review the relevant gamma-ray observations that can constrain the 
expected neutrino flux from the inner Galaxy. We suggest different possibili- 
ties to fit the gamma-ray data in scenarios where the emission is dominated 
by hadronic CR interactions. In Sec. [3] we use such fits to evaluate the ex- 
pected flux of neutrinos from the inner Galaxy. This has to be considered as 
the maximal neutrino flux which is compatible to all the available gamma-ray 
measurements. Any leptonic contribution to the observed gamma-ray emis- 
sion would lower the expectations. Our conclusions are provided in Sec. SJ 
The maximum level of the diffuse neutrino flux allowed by constraints from 
gamma-ray observations is within the detection capabilities of a km 3 -scale 
neutrino telescope located in the northern hemisphere. A detection would un- 
ambiguously reveal the hadronic origin of the multi-TeV emission detected by 
MILAGRO. 
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2 Diffuse gamma rays from the inner Galaxy 



The MILAGRO detection of a diffuse emission from both the Cygnus region 
(65° < I < 85°) and the more central region of the Galaxy (30° < I < 65°), rep- 
resents more than seven years worth of data collection. The large field of view 
of the detector (~2 sr) and its almost continuous data aquisition rate make 
the instrument highly suitable for discovering large-scale diffuse gamma-ray 
signals which Air-Cherenkov telescope instruments would find more difficult 
to detect. The most significant of the MILAGRO diffuse flux detections orig- 
inates from the Cygnus region and is 8.6 a above the background [T]. The 
central Galactic region was also detected with a smaller but still very robust 
significance of 5.1 a [2]. The median energy of the detected photons is in both 
cases 15 TeV. 

As here we are interested in the emission from the inner Galaxy, we plot in 
Fig. [1] the available gamma-ray measurement for two distinct regions of the 
Galactic plane. In the upper panel of Fig. [T] we show the MILAGRO point rel- 
ative for the region with Galactic coordinates 30° < / < 65° and —2° < b < 2°, 
together with the EGRET points for the same region. This is the region within 
the field of view of MILAGRO which is the closest to the Galactic centre. As 
mentioned in Sec. [1] and demonstrated in [2] , conventional models assuming 
that the gamma-ray emission is dominated by hadronic CR interactions, and 
that the CR spectrum has the canonical slope ~ 2.7 everywhere in the Galaxy 
fail to fit simultaneously the EGRET and MILAGRO data by roughly a factor 
of 5. A good fit to all the gamma-ray data can be obtained by assuming that 
the CR spectrum in the region considered here is not a single power law but 
rather a broken power law: 

m = Jo r E ~" s] e ( - E/Ecut) 



with slope a = 2 and a + S = 2.7 below and above E* = 80 GeV respectively. 
An exponential cutoff is added at an energy of E mt = 1 PeV, close to the 
energy of the CR knee. Jo is a free normalization parameter, determining 
the total energy in the form of CRs. A physical justification of such a spectral 
shape has been provided in [10J, where the authors suggested to interpret E* as 
the energy where the escape of CRs from the Galaxy changes from convective 
to diffusive. Thus, at low energies (below E*), where convection in the galactic 
wind dominates, the shape of the CR spectrum mimic the injection spectrum 
(« E~ 2 ), while at higher energies (above E*) diffusive escape of particle from 
the Galaxy steepens the spectrum to its canonical value ~ E~ 2 7 . 

The gamma-ray emission resulting from proton-proton interactions between 
CRs with the spectrum considered in Eq. [T] and interstellar gas has been 
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Fig. 1. Top panel: EGRET and MILAGRO data points for the region of the 
Galactic plane with coordinates 30° < I < 65°. The solid line represents a fit to 
both the sets of data due to 7r°-decay gamma-ray emission from the CR spectrum 
described by Eq. [TJ The dotted and dashed lines represents the gamma-ray emission 
from a hard (oc E~ 2 ) CR spectrum that fits the MILAGRO point only. Bottom 
panel: EGRET data points for the inner Galaxy (—40° < I < 40°). Lines are the 
same as in the top panel rescaled upwards of a factor ~ 1.3. 



calculated using the parameterization given in [20] and is plotted as a solid 
line in Fig. Q] (top paneljQ The curve fits both the EGRET and MILAGRO 
data, with the exception of the data points at energies below w 100 MeV, 
where electron Bremsstrahlung is well known to be the dominant mechanism 
for gamma-ray production in the Galactic disk pTJ|8] . 



A different possibility to explain the excess of multi-TeV diffuse emission above 



The best fit parameters slightly change if a different parameterization for the 
production of gamma rays is adopted (e.g. [21110122] ) since these parameterizations 
treat the production of low energy pions in proton-proton interactions in a different 
way. However, at the high energies which are most relevant for this paper all the 
approaches converge to the same spectrum, and the choice of the parameterization 
will not affect our final results. 
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the conventional models without invoking leptonic contributions is to assume 
the existence of a hard hadronic component in the gamma-ray spectrum that 
adds up to the canonical ~ E~ 21 component and dominates the total emission 
above « TeV energies. Such a hard gamma-ray spectra can be produced in 
different ways. One possibility is to have a number of unresolved CR sources 
within the emitting region, as suggested in p]. In this case, the CR spec- 
trum inside the sources is not steepened yet by diffusion in the interstellar 
medium and might thus contribute significantly to the multi-TeV emission. 
Hard gamma-ray spectra can also be produced by propagation effects when 
relatively young CR sources exist in an inhomogeneous interstellar medium 
which hosts molecular cloud complexes |24|I23| . This is a consequence of the 
energy-dependent nature of CR propagation away from their sources in the 
interstellar medium. Under certain conditions, an energy-dependent diffusion 
coefficient allows the formation of hard CR spectra around the sources, whose 
7r°-decay gamma-ray emission can be amplified by the presence of a massive 
molecular cloud. Thus, a population of CR accelerators (such as supernova 
remnants) might be responsible for the whole gamma-ray emission detected 
by MILAGRO from both the Cygnus and the inner Galactic region, which 
host dense molecular cloud complexes. 

In this paper we adopt a phenomenological approach, not focussing on any 
particular scenario for the production of the observed multi-TeV gamma rays. 
Being interested in obtaining the maximal flux of diffuse neutrinos we consider 
for further calculations two different CR spectra, which are felt to cover the set 
of plausible scenarios for a hadronic origin of the diffuse gamma-ray emission 
detected by MILAGRO: i) the broken power law described by Eq. [TJ which fits 
simultaneously the EGRET and MILAGRO data and ii) a hard CR spectrum 
of the form: 

J(E) = J E- 2 e { - E/E ^ , (2) 

which fits the MILAGRO point and is represented in Fig. [1] (top panel) by the 
dashed and dotted lines, characterized by two different values for the position 
of the high energy cutoff of 1 and 5 PeV respectively. The hard CR spectrum 
described by Eq. [2] provides a negligible contribution to the low energy (« 
GeV) gamma-ray emission detected by EGRET, that can be explained by 
conventional models that assume a steep CR spectrum [3ll8] . 

In order to estimate the diffuse neutrino flux from the inner Galaxy we pro- 
ceed as follows. We call here inner Galaxy the region of galactic coordinates 
-40° < / < 40° and -2° < b < 2°. This is the region for which the diffuse 
emission detected by EGRET shows the highest brightness (see Fig. 2a from 
[3]). Thus, if the emission has, as it is believed, a hadronic origin, the expected 
neutrino flux is also expected to be enhanced within this region. The missing 
piece of information is the gamma-ray flux from the inner Galaxy at multi-TeV 
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energies, which is the most relevant to constrain the neutrino emission. This 
flux is not available since the MILAGRO observations cover only a fraction 
of the inner Galaxy. We assume here that the enhancement detected by MI- 
LAGRO in the region of coordinates 30° < I < 65°, which partially overlaps 
with the region we consider, is representative for the whole inner Galaxy and 
that we can rescale it following the EGRET gamma-ray brightness profile. In 
Fig. [1] (bottom panel) we show the EGRET data points for the inner Galaxy, 
together with the three curves adopted in the upper panel to fit the gamma- 
ray data. All the curves have been multiplied by 1.3 in order to follow the 
spatial variation of the EGRET flux. This assumption is equivalent to saying 
that GeV and multi-TeV CRs are distributed in the inner Galaxy in the same 
way. It has to be kept in mind that the factor 1.3 takes into account of both 
the spatial distribution of CRs and the increase of the gas density toward the 
centre of the Galaxy. 

Before proceeding to estimate the expected neutrino flux, it is important to 
note that upcoming observations by GLAST will improve our knowledge of 
the Galactic diffuse gamma-ray emission and will extend the observed energy 
range up to ~ 100 GeV. Similarly, next generation instruments like HAWK 
will survey the gamma-ray sky in the energy range between ~ 1 and 100 TeV. 
The improved sensitivity and better coverage of the energy range will allow 
a more accurate fit to the data and consequently reduce the uncertainties in 
the estimate of the neutrino flux. 



3 Neutrino flux from the inner Galaxy 

In this section, after briefly reviewing the capabilities of km 3 -scale neutrino 
telescopes, we estimate the diffuse flux of neutrinos from the inner part of the 
Galaxy (-40° < / < 40°, -2° < b < 2°) by assuming that the CR spectra that 
we used in the previous section to fit the MILAGRO data are representative for 
the whole inner Galaxy. These are the CR spectra that produce the 7r°-decay 
emission plotted in the bottom panel of Fig. [IJ 

3.1 Capabilities of neutrino telescopes 

As a rule of thumb, km 3 -scale neutrino telescopes can detect a persistent and 
point-like source at a flux level of ~ 10 _11 z/ cm" 2 s" 1 after a few years of 
continuous observations. This flux is the total flux integrated above m 1 TeV 
and roughly corresponds to the observed flux of the Crab nebula in gamma 
rays. An accurate determination of the detection rate from a Crab-like source 
of neutrinos can be obtained by considering the telescope's effective area and 
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Fig. 2. Differential detection rates for neutrino sources with flux 
F(> 1 TeV) = W~ n u cm -2 s _1 . The dashed region represents atmospheric 
neutrinos. Left panel: Solid, dashed and dotted lines represent the count rates for 
neutrino sources with power law spectrum with slope 1.5, 2 and 2.5 respectively 
and high energy cutoff at E IPeV. Right panel: Count rates for neutrino 
sources with oc E~ 2 spectra and exponential cutoff at E cut = 10 and 100 TeV (solid 
and dashed line respectively). 

the spectrum of the neutrinos received by the source. The rate obtained in 
this way has to be compared with the detection rate of atmospheric neutrinos, 
which constitute the dominant background. All these aspects will determine 
which is the optimal energy range for the detection of astrophysical neutrinos. 

As an illustrative example, we consider here a hypothetical point-like and 
steady source of neutrinos with differential flux: 

with J normalised such that the integrated flux above 1 TeV is 10~ n cm' 2 s _1 . 
a and E cut are free parameters. For the effective area we used the one pro- 
vided in [25.26.27J for the KM3NeT detector. Being located in the northern 
hemisphere, this telescope will be able to observe the inner part of the Galaxy 
and thus is the more relevant for this paper. In the calculations, we used a 
convenient fit to the effective area which reads: 

A M E) = 10" E» ( °- 24 ) ( ° M ) (-*L-) cm* (4) 
V ; V0.24 + E7 \0.31 + EJ V37 + EJ V ; 

where E is the neutrino energy in TeV. The product between the effective area 
and the neutrino flux within one angular resolution element results in the ex- 
pected differential detection rate by KM3NeT, and it is shown in Fig. [2j Rates 
are defined as the number of neutrinos detected after one year of exposure. 

In the left panel of Fig. [2] we demonstrate the effect of changing the spectral 
slope of the neutrino source. The solid, dashed and dotted curves refer respec- 
tively to a = 1.5, 2 and 2.5, that cover the most plausible spectra for neutrino 
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sources. The position of the exponential cutoff is E cut ^> IPeV. In the right 
panel we fix a =2 and we change the energy at which the spectrum cuts off. 
Solid and dashed lines refer to E cut = 10 and 100 TeV respectively. In both 
the plots the shaded region represents the level of the atmospheric neutrino 
background, assumed to follow the spectrum predicted in [28]. The resulting 
atmospheric background rate is shown to be consistent with [29 30 J3T) . The 
adopted atmospheric background assumed a zenith angle of ~ 70°. This is 
the effective zenith angle of the inner Galactic region we are interested in 
(—40° < I < 40°, —2° < b < 2°), obtained by integrating the atmospheric 
neutrino flux over one day from source regions below the horizon. The energy 
resolution of the detector is assumed to be constant across the entire energy 
range. It is accounted for by convolving the contents of each energy bin with a 
Gaussian distributed in logarithmic energy space. The Gaussian has an RMS 
in the difference between the natural logarithm of the reconstructed and true 
neutrino energies of 0.3. The angular resolution of the detector is given by a 
parameterisation of the RMS difference between the reconstructed neutrino 
direction and the true neutrino direction from the same simulations used to 
produce the effective area curve [25|26] . and above 1 TeV drops below 0.2°. 
In these generic plots it is assumed the source is point-like, and both the sig- 
nal and background neutrino fluxes are integrated across a cone of optimal 
opening angle 1.58 times the angular resolution. 



It is evident from Fig. [2] that the prospects for km 3 size neutrino telescopes 
to detect such fluxes are optimal in the 10-1000 TeV energy range, where 
the signal flux is well above that due to atmospheric background events and 
the effective area grows sufficiently quickly with energy for detection to be 
possible. This result comes with the caveat that the cutoff energy of the flux, 
-E cut , sits at sufficiently high energies (>100 TeV). Thus the best targets for 
neutrino telescopes are sources exhibiting a hard spectrum which extends up 
to at least hundreds of TeV. It is for this reason that the detections of multi- 
TeV gamma-ray discrete sources and diffuse emission by MILAGRO are of 
such relevance for future prospects for the detection of the first high energy 
astrophysical neutrinos from the Galaxy. 



Qualitatively, similar conclusions can be reached also in the case of diffuse 
sources, the only difference being that we have to integrate both the signal 
and the atmospheric background over the extension of the sources. In the 
next section we will focus on this issue, by considering the expected neutrino 
emission from the inner part of the Galactic disk. 
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Fig. 3. Left panel: neutrino fluxes from the inner Galaxy. Solid, dashed and dot- 
ted lines refer to an underlying CR spectrum described by Eq. Q] and Eq. [2] with 
-Scut = 1 an d 5 PeV respectively. The flux of atmospheric neutrinos is shown as a 
shaded region. Right panel: Detection rates integrated above an energy E in a 
km 3 neutrino telescope. Curves have the same meaning as in the left panel. 

3.2 Neutrinos from the inner Galaxy: detectability for km 3 -scale detectors 



The left panel of Fig. [3] shows our predictions for the diffuse flux of neutri- 
nos expected from the inner Galaxy, defined by the coordinates —40° < I < 
40°, —2° < b < 2°. Neutrino spectra are calculated following [20J and are 
the sum of all the neutrino flavours produced in proton-proton interactions 
(two muon neutrinos plus one electron neutrino) divided by three to take into 
account neutrino oscillations (maximal mixing). The adopted parent proton 
spectra are the same we used to fit the MILAGRO data (see Sec. [2]) and then 
extrapolated to the whole inner Galaxy, namely, we used the same CR proton 
spectra that produce the gamma-ray emission plotted in the bottom panel of 

Fig. m 

The solid line represents the neutrinos produced by CRs with spectrum de- 
scribed by the broken power law in Eq. [TJ while dashed and dotted lines refer to 
a hard CR spectrum like in Eq. [2j with exponential cutoff at energies E cut = 1 
and 5 PeV respectively. The shaded region represents the atmospheric neutri- 
nos, integrated over the whole of the considered region, which has an extension 
of 0.097 steradians. 

For the case of the broken power law the atmospheric neutrinos dominates over 
the signal for all neutrino energies. On the other hand, if a hard spectrum is 
considered, and the exponential cutoff in the CR spectrum is significantly 
above 1 PeV, the neutrinos from CR interactions in the inner Galaxy start to 
dominate over the atmospheric background above an energy of a few tens of 
TeV. 

In order to investigate the capability of a neutrino telescope to detect the 
diffuse emission we convolved the neutrino fluxes with the effective area (Eq. HJ) 
to obtain the integral detection rate. The rates above a given energy E for one 
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year of exposure are shown in the right panel of Fig. [3j We neglect here the 
effect of a finite energy resolution of the detector. Line types are the same as in 
the left panel of Fig. [3J It is evident that the detection of diffuse neutrinos from 
the inner Galaxy becomes most feasible for neutrino energies above a few tens 
of TeV, where the signal becomes comparable or above the background counts 
rate. Above 10 TeV a rate of 5, 9 and 15 neutrinos is expected after a single 
year of exposure for the solid, dashed and dotted curve respectively, against 
a background of about 28 neutrinos per year. For the most optimistic case 
(dotted line in Fig. [3]) the neutrino signal becomes equal to the background 
at an energy of 20 TeV with about 10 events per year. 

The visibility of the inner Galaxy must also be considered, which for a detec- 
tor located in the Mediteranean at 42°50W, 6°10'£, is visible for « 70% of 
the year. The rates shown in Fig. [3] have to be considered as the number of 
neutrinos after one year of exposure, which translates in w 1.5 yr of physical 
time. 



4 Discussion and conclusions 

In this paper we presented an estimate of the neutrino flux from the inner 
Galaxy, defined as the region of coordinates —40° < I < 40° and — 2° < b < 2°. 
Our calculations rely on the following assumptions: 

(1) the diffuse gamma-ray emission detected at multi-TeV energies by MILA- 
GRO from the region of the Galactic plane of coordinates 30° < I < 65° 
is representative for the whole inner Galaxy; 

(2) the diffuse galactic GeV and multi-TeV emission detected by EGRET 
and MILAGRO respectively follows the same spatial distribution; 

(3) the diffuse gamma-ray emission is entirely of hadronic origin. 

These assumptions allow to obtain an almost model independent estimate for 
the diffuse flux of neutrinos from the inner Galaxy. Such an estimate has to be 
considered as the maximal expected flux of neutrinos compatible with gamma 
ray observations, since any leptonic contribution to the observed gamma ray 
emission would lower our expectations. 

The maximal expected neutrino flux is within the detection capabilities of 
a km 3 -scale neutrino telescope located in the northern hemisphere. For the 
most favorable assumptions made in Sec. [3J we expect from the inner Galactic 
region 15 neutrinos above 10 TeV per year of exposure against a background 
of 28 atmospheric neutrinos. At 20 TeV the signal equates the background 
with about 10 events per year. The detection of these neutrinos would consti- 
tute an unambiguous proof of the hadronic origin of the gamma-ray emission 
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detected by MILAGRO. Finally, it is important to stress that MILAGRO is 
operating in an energy range similar to the one where neutrino telescopes reach 
their best performances. This makes MILAGRO and future instruments like 
HAWK, the ideal gamma-ray detectors for multi-messenger studies of sources 
of astrophysical neutrinos. 
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